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Evaporating spray plays an important role in spray combustion processes. This paper describes the 
development of a new finite-conductivity evaporation model, based on the two-temperature film theory, for 
two-phase numerical simulation using Eulerian-Lagrangian method. The model is a natural extension of the 
T-blob/T-TAB atomization/spray model 1 which supplies the turbulence characteristics for estimating 
effective thermal diffusivity within the droplet phase. Both one-way and two-way coupled calculations were 
performed to investigate the performance of this model. Validation results indicate the superiority of the 
finite-conductivity model in low speed parallel flow evaporating sprays. High speed cross flow spray results 
indicate the effectiveness of the T-blob/T-TAB model and point to the needed improvements in high speed 
evaporating spray modeling. 


Nomenclature 


A d ■= Droplet surface area (m 2 ) 

B m = Gas phase mass transfer number (Spalding number) 

C d = Discharge coefficient of injection nozzle (= 0.7) 

C p = Specific heat capacity of mix (KJ/kg.K) 

Cp - Turbulence constant 

D = Diffusivity (m 2 /s) 

Kc = Loss coefficient due to nozzle inlet geometry (= 0.45) 

K e = Proportionality constant (= 0.23) 

k\ - Turbulent kinetic energy of liquid (m 2 /s 2 ) 

L v = Latent heat of the fuel at the surface temperature (KJ/kg) 

m* = Mass of droplet (kg) 

rh d = Droplet mass evaporation rate (kg/s) 

Nu g = Gas phase Nusselt number 

Pr t = Turbulent Prandtl number (=0.9) 

q g = Heat transfer rate from gas phase to droplet surface 

qi = Heat transfer rate from droplet surface to droplet core 

R p = Droplet radius (m) 

Sh g = Gas phase Sherwood number 

s = Area ratio at nozzle contraction (= 0.0 for this study) 

T g = Gas temperature of environment (K) 

T d = Bulk temperature of droplet (K) 
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T s = Surface temperature of droplet (K) 
t = Time (s) 

ai am = Laminar thermal diffusivity (m 2 /s) 

ctturb = Turbulent thermal diffusivity (m 2 /s) 
ctgff = Effective thermal diffusivity (m 2 /s) 

Xj = Liquid thermal conductivity (KW/m.K) 

5 e = Equivalent thickness of thermal boundary layer (m) 

E\ = Turbulent dissipation rate of kinetic energy of liquid (m 2 /s 3 ) 

p = Density (kg/m 3 ) 


L INTRODUCTION 

I n liquid fueled combustion systems such as rocket engines and ramjet/scramjet engine applications, the 
combustion performance depends strongly on liquid atomization, spray vaporization, and the mixing process 
between the free stream environment and the liquid fuel spray. Smaller drop sizes generated from the spray devices 
increase the specific heat and mass transfer area of the fuel and thereby achieve high rates of mixing and 
evaporation. Consequently, higher engine efficiency would be realized in the combustion process. A large diversity 
of multi-phase gas-liquid flows of practical interest involves the evaporation of near spherical liquid droplets in high 
temperature turbulent environments. Such flows cover a wide range of applications including spray combustion. 
These situations involve a dispersed liquid phase species in the form of a large number of discrete droplets 
converting and vaporizing in a continuous gas phase species. As droplets enter the hot free stream, only a small 
fraction of the liquid droplet near the liquid/gas interface heats up while the core region remains cold. The heat 
transfer through conduction and convection (possibly by internal turbulence and circulation) to its interior proceeds 
until the end of its lifetime. Resolving the transient phenomena within liquid droplets, as well as the interfacial 
heat/mass transfer are the subjects of various modeling approaches 3-13 . These evaporation models, when 
implemented in modem computational analyses utilizing the Eulerian-Lagrangian method, would involve 
calculations of velocity, temperature for both and vapor of the surrounding gases and the liquid droplets, as well as 
species mass fractions within the gas phase. 

Most of the current droplet evaporation models used for computational analysis of spray combustion 
phenomena is based on the classical constant droplet-temperature model (D 2 law) or infinite conductivity (perfect 
mixing within the liquid droplets) models. Even the performance of the gas-side vaporization models may vary, 
droplet evaporation rates predicted by liquid-side infinite conductivity models tends to over predict or under predict 
the evaporation mass flux, depending on the ambient temperature conditions. Bertoli and Migliaccio 3 showed that 
the accuracy of CFD computations of heating, evaporation and combustion of diesel fuel sprays could be 
substantially increased if the assumption of infinitely high thermal conductivity of liquid is relaxed. Sprays exposed 
to convective flows behave very differently from quiescent clouds of drops. Several modeling approaches have been 
proposed to account for finite conductivity effects within the evaporating droplets. Based on single drop analysis, 
Sirignano and Abramzon 4-5 takes into account details of external flow around a droplet and the coupled internal flow 
dynamics inside the droplet for a droplet immersed in a convective flow. However these results are not applicable 
for dense spray conditions. Sazhin et al. 6-7 pointed out that it is not economical (for CFD spray combustion) to 
resolve temperature profiles within evaporating droplets and proposed an alternative to numerically modeling the 
fuel droplet heating by obtaining analytical solutions to the heat conduction equation within the droplet. Miller et al. 9 
introduced non-equilibrium Langmuir-Knudsen evaporation law for use in droplet combustion models and found 
non-equilibrium effects to be important for droplet sizes found in practical spray calculations. Renksizbulut et al. 10 
developed a Nusselt number vaporization model which incorporated heat transfer enhancement due to internal 
circulation. Reveillion and Vervisch 10 have considered the effects of both internal circulation and mixing inside the 
droplet. Takahashi et al. 11 modeled a turbulent thermal conductivity for the heat transfer behavior in the sheet region 
based on the model developed by Celata et al. 12 

The goal of this study is to extend a recently developed atomization/spray model 1-2 , the T-blob/T-TAB model, 
to include spray evaporation effects. Due to the unique feature of T-blob/T-TAB, in which the turbulence 
characteristics is accounted for within the droplet phase, extension of this model to include finite conductivity effect 
in the evaporating droplet can be made naturally. The model development will be described in this paper. 
Preliminary validations for one-way coupling, two-way coupling 2-D axis-symmetric and 3-D evaporating sprays 
will be presented with a long term goal of developing numerical capabilities for both low speed and high speed 
spray combustion applications including liquid turbulence effects. 
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II. NUMERICAL MODELING APPROACHES 


2.1 T-blob and T-TAB Atomization/Spray Models 

The current vaporization model is developed for computational analysis based on the Eulerian-Lagrangian 
numerical approach. In this formulation, the spray/dropiets dynamics is described in a Lagrangian coordinate such 
that numerical droplets are tracked within the Eulerian gas dynamics. Liquid phase is tracked from the injector 
plane, and the primary atomization, as well as the subsequent secondary breakup is modeled using the T-blob/T- 
TAB hybrid model. Both primary and secondary droplet breakup processes are modeled and the transition from 
primary to secondary breakup is modeled based on energy balance. In addition to the droplets position and velocity, 
liquid turbulence is accounted for through the injector characteristics by the two-equation k-e turbulence model 
formulation using the T-blob/T-TAB model. Detailed model description and validations can be found in references 1 
and 2, and it is suffice to say that within each numerical droplet, turbulence characteristics such as fluctuating 
velocity level, length and time scales are supplied by the model. 

2.2 The Finite Conductivity Evaporation Model 

To relax perfect mixing assumption within the liquid droplet, a “two-temperature model” is formulated. This 
type of two-temperature modeling has been described in the non-equilibrium evaporation models by Miller et al. 8 , 
and was used by Ra and Reitz 13 for gasoline engine applications. In the two-temperature model, the core (or bulk) 
temperature (T d ), is assumed well-mixed by convection/turbulence transport. In consistence with the “film theory”, 
heat resistance exists at the near surface region, and the droplet surface temperature (T s ) differs from the droplet core 
temperature. The heat transfer coefficient across this thin film (or boundary layer) is then formulated through the 
turbulence characteristics supplied from the T-blob/T-TAB model, to account for the finite conductivity effect. In 
the Lagrangian coordinate, the heat-up of the droplet core is formulated as: 


where h t is the liquid-side heat transfer coefficient and A d is the droplet surface area. The heat transfer coefficient 
is determined from an effective thermal diffiisivity formulation 5 as; 





( 2 ) 


where X 7 is the liquid thermal conductivity and the film thickness, S e , is 


given 


by 


eff t . The effective 


thermal diffiisivity ( a eff ), based on the turbulence characteristics within the droplet, is estimated from: 
a ejf =a !am + a mrh , in which the turbulent thermal diffiisivity is calculated from the two-equation turbulence model 


diffiisivity formulation: a 


turb 


Pr, 


, here Pr f is the turbulent Prandtl number. The liquid droplet turbulence 


quantities k { and £ ( are obtained from the T-blob/T-TAB atomization/spray model. 

In cases where gas phase diffusivities are much larger than liquid diffusivities, the droplet core heating will be 
rate controlling and the gas-side heat/mass transfer will respond in a quasi-steady manner. The surface temperature 
of the droplet is determined from a heat and mass transfer balance at the interface between the droplet and the 
surrounding gas assuming no heat accumulation at the droplet surface such that: 


L v 


m d =q g ~q,. 


( 3 ) 


where L v - latent heat of the fuel at the surface temperature (KJ/kg), 


q g - heat transfer rate from the environmental gas to the surface. In this paper, the classical Spalding evaporation 
model is used to model the gas-phase transport, thus the gas heat transfer rate was calculated as: 
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= 2 xD p Z g Nu g (T g -Z)!2^ 

2 (T — T 1 

q J - heat transfer rate from the droplet interior to droplet surface = JX-± a , and the evaporation rate at 

5 e 

the surface is given as: 

= 2*D p (p g D)Sh g ln(l + B m ) (4) 

In the above equation, B m is the Spalding mass transfer number, Sh g is the Sherwood number, and D is the binary 
diffusivity. The Sherwood and the Nusselt numbers were calculated using the classical correlations, given by: 

Nu g = 1 + 0.3 Re p ° 5 Pr 0333 (5) 

Sh g =1 + 0.3 Re p ° 5 Sc 0 333 (6) 

The equilibrium condition at the droplet surface is determined by the Clausius-Clayperon relationship. 


III. VALIDATION RESULTS AND DISCUSSION 


Several test cases were selected to investigate the performance of the current two-temperature evaporation 
model. Initially, a simple FORTRAN program was written for one-way coupled evaporating atomizing spray in a 
quiescent gas. This program is an extension of the T-blob/T-TAB code written by Trinh and Chen 1 . Detailed 
numerical calculation issues such as time step, initial conditions and the iterative loop of obtaining droplet surface 
temperature and core temperatures can be found in Balasubramanyam 14 . After the initial testing, the current model 
then was implemented into an existing CFD code, CFD-ACE+ 15 , for two-way coupling Eulerian-Lagrangian 
calculations. Both 2-D axis-symmetric co-flowing cases and 3-D cross flow cases were investigated. Again, Ref. 14 
should be consulted for detailed numerical (such as grids and time steps) implementation issues. 


3.1 One-Way Coupled Evaporating Atomizing Spray 


This case is an extension of the one-way T-blob/T-TAB testing case described in Trinh and Chen 1 . The purpose 
here is to investigate the concept of boundary layer film thickness within with liquid droplets involving two 
temperatures. In the course of study, it was found that, due to the fact, in the secondary breakup regime (i.e. the T- 
TAB regime) the droplets were so small they were heated up rapidly. Thus the current model was only implemented 
within the T-blob (i.e., the primary breakup) model. The Tridecane fuel was issued through a long injector tube at 
300°K. The length of the injector nozzle (L) is 0.8 mm and the nozzle diameter ( d nozzle ) is 0.3 mm. A jet injection 


velocity (U) of 102 m/s was used for one of the test cases. Based on these values, the initial conditions for liquid 
turbulent kinetic energy and its dissipation rate are estimated 1 as: 

U! 


81 /D, 

, K — 
Ke 2L 


nozzle 


(7) 


i-K.-O-s 1 ) 


( 8 ) 


The environment is quiescent nitrogen at a temperature of 600°K. The gas properties were calculated based on 
the reference state determined from the ‘ l/3 rd ’ rule 5 . In this calculation, a ‘blob’ of numerical droplet was injected at 
the orifice plane with orifice diameter. The droplet then went through first and secondary breakup processes, thus its 
diameter decreased in time. The variation of the thermal boundary layer within the droplet also changed in time. In 
figure 1, the change of the thermal boundary layer normalized with the droplet radius versus time is shown for the 
test case of a droplet evaporating in a quiescent environment. It was also observed that the turbulent diffusivity 
evaluated using the dispersed phase flow turbulence was at least two orders of magnitude higher than the laminar 
thermal diffusivity. The result shows that the normalized thermal boundary layer is rather small and exhibits 
reasonable physical trend. Detailed droplet characteristic history, such as the two temperature history, heat up time, 
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as well as turbulence characteristics can be found in Ref. 14. The model was then incorporated into the finite volume 
commercial CFD code ACE+ 15 . Current work is focused towards validation of the code with test cases for which 
experimentally measured data is available for comparison. Some of these results will be presented in this paper. 

3.2 Two-Way Coupled Evaporating Spray Validation 

To evaluate the current evaporation model, the 2-D axis-symmetric subsonic low-speed evaporating spray of 
Yakota et al. 16 was tested first for a two-way coupling CFD calculation since a similar non-evaporating was used for 
the T-blob/T-TAB validation study 2 . Liquid fuel (Tridecane Ci 3 H 28 ) is injected through a single-hole nozzle into a 
high pressure, high temperature ambient N 2 environment. The initial test conditions for the evaporating spray are 
summarized in Table 1 . 


Table 1. Test Conditions for the Measurement of Yokota et al. 16 


Case 

mm 

Pgas 

T a mb 

Mil >i 

Gas Environment 

Y-l 

30 Mpa 

3 Mpa 

900 K 

0.00326 kg/s 

n 2 


The nozzle diameter was 0.16 mm. Equations (7) and (8), with the same values for nozzle parameters except for 
nozzle diameter and jet velocity, were used to estimate the initial liquid turbulence quantities. A computational 
domain of 20 mm in radius and 100 mm in length discretized by a 50 radial x 75 axial grid was used. The mesh 
spacing was non-uniform with refinement on the centerline and close to the injector. A constant time step of 2.5E-6 
sec. was used with an injection period of 4 psec. The properties of liquid fuel Tridecane were taken from the 
NIST/JANAF database. Estimating the penetration of a fuel jet into an air stream is an important global property for 
model validation and is presented first. In Figure 2, the predicted tip penetration results using the current finite- 
conductivity (F-C) evaporation model coupled with the T-blob/T-TAB atomization model are compared with the 
measured data. For reference, predictions using classical Blob/TAB/infinite-conductivity (I-C) model, as well as 
using T-blob/T-TAB/infinite-conductivity model are also shown in the same figure. It can be seen that the infinite- 
conductivity evaporation models tend to over predict the evaporation rate, thus gave shorter tip penetration. On the 
other hand, the finite-conductivity model slows down the evaporation process, and produce larger droplets and 
longer penetration. In figure 3.b, it can be observed that larger numerical droplets are clustered near the tip region 
due to slower evaporation process calculated by the finite-conductivity model, when compared with the parcel 
distribution using infinite-conductivity model 3. a, 1 psec after the injection. It should be noted that the coalescence 
model 2 was used for all two-way coupling calculation cases. The coalescence model is responsible for the calculated 
overshoot phenomena observed in the initial period of injection for all simulated cases. 

The heat and mass transfer aspects of the evaporating jet are shown in figures 4 and 5 at time of 4 psec after 
injection. As can be observed, the models incorporating the turbulence effects in the primary and secondary 
atomization processes give more reasonable qualitative pictures when compared with the classical atomization 
model without liquid turbulence effect. Comparing figures 4.b and 4.c, the surface temperature contours predicted 
by the current finite-conductivity model show the effect of slowing down the rate of evaporation, based on the 
turbulence levels experienced by the individual droplets, rather than the constant drop temperature assumption of the 
classical infinite conductivity model. The fuel mass fraction contours for test case Y-l are as indicated in figures 5.a, 
5.b and 5.c. The comparisons between figures 5.b and 5.c show that diffusion of the fuel species into the 
environment is at much lower rate than the case with the classical infinite conductivity model. 

3.3 Two-Way Coupled 3-D Supersonic Cross-Flow Validation 

The study of liquid atomization process in a high-speed flow has become an important area of research in air- 
breathing vehicle (Scramjet) development. The breakup process in a cross-flow involve the effects of drop 
deformation, displacement, breakup time and ignition were studied by Ranger and Nichols 19 and Shetz et al. 26 . 
Though advances (Gruenig and Mayinger 21 , Yu et al. 22 , Vinogradov et al. 23 , Inamura et al. 17 ) in this area of research 
have been made it is still not progressed so far as to provide validation cases for evaporation measurements due to 
the complexity of the phenomena occurring in these types of situations and the lack of measurement techniques to 
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capture these phenomena. Owing to limitations in availability of data for comparisons it was decided to apply the T- 
blob/T-TAB (originally developed for low speed parallel injector flows) model coupled with the infinite 
conductivity model, for simulating the test case of Lin et al. 24 (also used by Im et al. 18 ) for the injection of a water jet 
into a cross-flow at Mach 1 .94. 

A three dimensional computational domain with dimensions of 762x152x127 mm, discretized with a 90x25x30 
grid distribution, with refinement at the injection point was used for the numerical computations. The inlet boundary 
conditions for the chamber were given by a static pressure of 29 KPa, a static temperature of 304. 1°K and a flow 
velocity of 679 m/s based on a free stream Mach number of 1.94. The injection velocity of the spray jet was 
determined based on a jet-to-air momentum ratio q 0 = 7.0. In this study the water jet injection temperature was set to 
300°K. The injection nozzle diameter was 0.5 mm and located at a distance of 139 mm from the chamber inlet. 
Again, inlet turbulence quantities are estimated by Eqs. (7) and (8) due to the lack of information. A computational 
time step of 2.5E-6 jxsec was used with an injection period of 4 psec. 

The cross-section of the jets is circular at the injector exit and becomes kidney shaped due to the dynamic 
pressure of the air stream. This phenomenon was observed by Shetz and Thomas 25 . The spray tip penetration for this 
cross-flow test case is shown in figure 6. The jet penetration increased rapidly in the vicinity of the injector exit and 
then gradually increased due to increase in the drag of the air stream. The momentum of the spray droplets increases 
due to acceleration by the airflow in the direction of the supersonic airflow. Consequently, the penetration curve of 
the spray plume approaches the airflow direction as observed in figure 7. Figures 8 and 9 show the Mach and 
pressure contours at the mid-plane cut of the evaporating spray. An experimental correlation curve suggested by Lin 
et .al, 24 is also shown in figure 6, It should be pointed out that the current simulation employs the same T-blob/T- 
TAB model and the associated low-speed droplet drag law without adjusting any model constants. However, the 
droplets have been observed to accelerate to supersonic speeds at a later part of the spray. The evaporation of 
droplets would also result in a lowering of the drag force of the air stream 12 . Further improvement in drag law to 
account for high speed effect is definitely needed. In addition, the inlets liquid jet turbulence level was found to be 
quite sensitive in terms of cross-flow spray penetration calculation 14 . Recently, inlet liquid jet turbulence is also 
observed to be sensitive to spray dynamics experimentally 21 . Considering the uncertainties associated with the jet 
inlet conditions, the comparison with experimental data is reasonable. Further numerical studies should be 
conducted to clarify the effects of air temperature on jet penetration incorporating the finite-conductivity 
evaporation model. 


IV, CONCLUDING REMARKS 

We have developed a new finite-conductivity model, based on the two-temperature film theory, for evaporating 
spray numerical calculations. The model is a natural extension of the T-blob/T-TAB atomization/spray model which 
supplies the turbulence characteristics for estimating effective thermal diffusivity within the droplet. Based on one- 
way coupling simple spray results, the model exhibits reasonable physical trends in terms of droplets evaporation 
features. In two-way Eulerian-Lagrangian multi-dimensional full CFD simulations, the current finite-conductivity 
model coupled with the T-blob/T-TAB model show superior performance to conventional infinite-conductivity 
evaporation models when comparing to low-speed parallel injector flow experimental data. Due to the lack of high 
speed evaporating spray data, limited results using the same T-blob/T-TAB model for high speed cross flow spray 
were presented. Extension of the low-speed T-blob/T-TAB model for high-speed cross flow application show 
surprisingly reasonable results considering the uncertainties involved. The preliminary high speed evaporating spray 
study also pointed to the future directions of improving the current model in terms of high speed droplet drag laws, 
additional breakup mechanisms, and needs of detailed quality experimental data. 
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Figure 1. Variation of the normalized thermal boundary layer with time. 





Figure 3. Spray Parcel Distribution with T-blob and TTAB at 1 psec; (a) I-C Model 

(b) F-C Model 
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Figure 4. Temperature Contours at 4 psec; (a) I-C with Blob/TAB Model, 
(b) I-C with T-blob /TTAB, (c) F-C with T-blob /TTAB 
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Figure 5. Fuel Mass Fraction Contours at 4 jisec; (a) I-C with Blob/TAB Model, 
(b) I-C with T-blob /TTAB, (c) F-C with T-blob/TTAB Model 
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Figure 7. Jet penetration pattern at mid plane 
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Figure 8. Mach number contours - I-C with T-Blob/T-TAB Model 



Figure 9. Pressure contours - I-C with T-Blob/T-TAB Model 
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